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Abstract 

-Eg X Eg heterotic string and M-theory, when appropriately com- 
pactified, can give rise to reahstic, = 1 supersymmetric particle 
physics. In particular, the exact matter spectrum of the MSSM, 
including three right-handed neutrino supermultiplets, one per fam- 
ily, and one pair of Higgs-Higgs conjugate superfields is obtained by 
compactifying on Calabi-Yau manifolds admitting specific S'[/(4) vec- 
tor bundles. These "heterotic standard models" have the SU{3)c x 
SU{2)l X U{1)y gauge group of the standard model augmented by 
an additional gauged U{1)b-l- Their minimal content requires that 
the B-L gauge symmetry be spontaneously broken by a vacuum ex- 
pectation value of at least one right-handed sneutrino. In a previous 
paper, we presented the results of a renormalization group analysis 
showing that B-L gauge symmetry is indeed radiatively broken with 
a S-L/electroweak hierarchy of 0(10) to 0(10^). In this paper, we 
present the details of that analysis, extending the results to include 
higher order terms in tan(3~^ and the explicit spectrum of all squarks 
and sleptons. 



1 Introduction 



Smooth compactifications of the weakly coupled [T] and strongly coupled 
[21 [3l H] i^s X i?8 heterotic string have been studied for many years. When 
the compactification is on a Calabi-Yau threefold with a slope-stable, holo- 
morphic vector bundle, the low energy four-dimensional effective theory is 
N = 1 supersymmetric. In recent years, such compactifications have been 
extended to complete intersection and elliptically fibered Calabi-Yau spaces 
admitting vector bundles constructed using monads [3 [6l [71 El [9], spectral 
covers [TOl [TTl [121 [13] and by extension of lower rank bundles [HI [15] . The 
formalism for explicitly computing the low energy spectrum in each of these 
cases has been developed, and presented in [THl [H], [ISl [IS] and [201 [21] 
respectively. Cohomological methods have been used to calculate the tex- 
ture of Yukawa couplings and other parameters in these contexts [221 [231 [21] • 
Finally, the non-perturbative string instanton contributions to the super- 
potential have been computed [2S1 [2S1 [23 [2H] and used to discuss moduli 
stability, supersymmetry breaking and the cosmological constant [22] • These 
methods underlie the theory of "brane universes" [301 [31] and new approaches 
to cosmology [321 [33] • 

In a series of papers, compactifications of the Eg x Eg superstring have 
been constructed on elliptically fibered Calabi-Yau spaces with Z3 x Z3 ho- 
motopy over a dPg surface |31]. These spaces admit a specific class of slope- 
stable, holomorphic vector bundles with structure group SU (4) that are con- 
structed by extension and are equivariant under Z3 x Z3 [35] . The non-trivial 
homotopy allows one to extend these bundles with flat Z3 x Z3 Wilson lines. 
Using the methods referenced above, the four- dimensional spectrum of these 
theories was computed. It is found to have precisely the matter content of 
the minimal supersymmetric standard model (MSSM), including three right- 
handed neutrino chiral supermultiplets, one per family. In addition there are 
a small number of Higgs-Higgs conjugate pairs, two in the model presented 
in [36] and one in the vacuum of [37]. These are termed "heterotic standard 
models" . They all contain a relatively small number of geometric and vector 
bundle moduli and each possesses an acceptable texture of Yukawa couplings. 
For specificity, in this paper we will consider the minimal heterotic standard 
model [37] containing the exact matter spectrum of the MSSM, one pair of 
Higgs-Higgs conjugate superfields as well as three complex structure moduli, 
three Kahler moduli and thirteen vector bundle moduli. 

The four-dimensional gauge group is obtained through the sequential 
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breaking of Eg by the S'f/(4) structure group of the vector bundle and the 
Z3 X Z3 of the Wilson lines. We find that 

Es '-^^ Sptn{10) SU{3)c x SU{2)l x U{l)y x U{1)b-l • 

Note that in addition to the standard model gauge group SU (3)c x SU (2) l x 
U{1)y, there is an extra gauged U{1)b-l- This arises from the fact that 
Spin{10) has rank five and that the rank must be preserved when the group 
is further broken by any Abelian finite group, such as Z3 x Z3. Since the 
standard model group has rank four, an extra U{1) gauge factor must appear, 
in this case precisely U{1)b^l- 

We would like to emphasize that there is a direct relationship between 
having three right-handed neutrino chiral multiplets, one per family, and the 
appearance of the additional U{1)b-l gauge factor. The S'f/(4) structure 
group in heterotic standard models is chosen precisely because the decompo- 
sition of the 248 of Eg with respect to it contains the 16 representation of 
Spin{10). It is well-known that each 16 is composed of one family of quarks 
and leptons, including a right-handed neutrino. This fact makes choosing a 
vector bundle with an SU (4) structure group a natural way to ensure that 
the spectrum contains three right-handed neutrinos. However, the rank of 
Spin{10) is five, one larger than the standard model gauge group. Hence, 
when Spin{10) is broken to SU{3)c x SU{2) l x U{1)y by the Abelian Z3 x Z3 
Wilson lines, an additional U{1)b-l must appear. 

Interestingly, the appearance of this U{1)b-l factor naturally solves an- 
other important phenomenological problem. As is well-known, the MSSM 
must be augmented by a specific Z2 symmetry called i?-parity (or matter par- 
ity, its supersymmetric extension) [381 EHl HOl SI] • This prohibits dangerous 
baryon and lepton number violating operators that would lead, for example, 
to rapid nucleon decay. Unfortunately, it has proven difficult to obtain dis- 
crete matter parity as a symmetry of smooth heterotic compactifications. As 
we have seen, however, a gauged U{1)b~l symmetry arises naturally in com- 
pactifications with 5*^7(4) structure group. U{1)b^l contains matter parity 
as a subgroup and, hence, as long as U{1)b-l is unbroken, matter parity is a 
symmetry. Since U{1)b-l is not observed in electroweak processes, it must 
be spontaneously broken at a higher scale. In supersymmetric grand unified 
theories (GUTs), where the matter content can be chosen arbitrarily, this 
breaking can occur at a high scale [121 HHl HH |16]. This is accomplished 
by adding multiplets to the MSSM, neutral under the standard model group. 
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for which 3{B — L) is an even, non-zero integer [17]. It is then arranged for at 
least one of these multiplets to get a large vacuum expectation value (VEV). 
This spontaneously breaks the gauged symmetry, giving a large mass to the 
B-L vector boson, but, since the 3{B — L) charge is even, leaves Z2 matter 
parity as a discrete symmetry. Hence, although U{1)b~l is broken at a high 
scale, Z2 matter parity is an exact symmetry at lower scales. However, this 
mechanism cannot occur within the context of the exact MSSM spectrum in 
the heterotic standard models, since for all fields 3{B — L) is either ±1, ±3 
or 0. It follows that when U{1)b-l is broken by a non-zero VEV, there is 
no residual Z2 symmetry. Therefore, in heterotic standard models U{1)b-l 
gauge symmetry must be broken at a low scale, an order of magnitude or 
two above the electroweak scale. But how can this be accomplished? 

In the MSSM spectrum, the only scalar fields that carry a non-trivial 
B-L charge, but transform trivally under the standard model gauge group, 
are the right-handed sneutrinos. Thus, B-L symmetry must be broken by 
at least one of these scalars acquiring a non-vanishing VEV from radiative 
corrections. However, to sufficiently suppress large baryon and lepton number 
violation, this must occur just above the electroweak scale. To analyze this, 
one adds to the supersymmetric MSSM the soft supersymmetry violating 
operators that arise from various sources, such as gaugino condensation and 
the moduli vacuum state, during compactification. Whatever the source, 
these operators are of a specific form first worked out in [IS] and discussed 
within the context of generalized scenarios in [5U| \5T[ [33] . The initial 
values of the parameters are set by the details of the compactification, and 
are generically moduli dependent. At any lower scale, these parameters are 
determined by a complicated set of intertwined, non-linear renormalization 
group equations (RGEs) [Ml ESI ESI EZl ESI EHl ED]- It is by no means clear 
that a non- vanishing sneutrino VEV will necessarily develop. If it does, one 
must still show that a neutral Higgs field will get a non-zero VEV, thus 
breaking electroweak symmetry, at a scale an order of magnitude or two 
lower than the sneutrino VEV. Finally, it is of interest to know whether this 
result requires extremely fine-tuned parameters or is, more or less, a natural 
hierarchy. 

In a previous paper [61|, we presented the results of a renormalization 
group analysis of the minimal heterotic standard model with a reasonable 
set of assumptions about the initial soft supersymmetry breaking param- 
eters. These assumptions are consistent with basic requirements of phe- 
nomenology, such as suppressed fiavor changing neutral currents, but are 



3 



further constrained so as to allow a quasi-analytic solution of the RGEs. We 
found that B-L symmetry is indeed spontaneously broken by a radiatively 
induced VEV of at least one right-handed sneutrino. Electroweak symmetry 
is then radiatively broken by a Higgs VEV at a lower scale, with the B- 
L/electroweak hierarchy of 0(10) to O(IO^). The purpose of this paper is to 
present the detailed renormalization group calculations leading to those con- 
clusions. These include both analytic, quasi-analytic and purely numerical 
solutions of the relevant equations. 

Specifically, we will do the following. In Section 2, the chiral fields of the 
U{1)b-l extended MSSM are presented and their supersymmetric interac- 
tions, via the superpotential and D-terms, are discussed. The complete set of 
soft supersymmetry breaking operators in this context are then introduced. 
Section 3 is devoted to presenting and solving the RGEs associated with the 
spontaneous breaking of the U{1) b-l gauge symmetry. This analysis involves 
the gauge parameters, gaugino and slepton masses and both the B-L and Y 
Fayet-Iliopoulos (FI) parameters [62]. Using these results, it is shown that 
U{1)b-l is indeed radiatively broken by a non-zero right-handed sneutrino 
VEV, and the details of this vacuum are presented. In Section 4, the analysis 
is extended to include both up and down Higgs and squark masses, as well 
as the /i and B parameters. It is then shown that at the B-L scale, elec- 
troweak symmetry, as well as color and charge, remain unbroken. All RGEs 
are then scaled down several orders of magnitude. We demonstrate that a 
Higgs VEV now develops which spontaneously breaks electroweak symmetry, 
without breaking color or charge. The B-L and Higgs VEVs are presented 
and all squark, slepton and Higgs masses are calculated in this vacuum. The 
results are a detailed function of the initial right-handed sneutrino and Higgs 
mass parameters, m,^(0) and mH{Q) respectively, as well as inverse powers of 
tanfj. The relationship between mj,(0) and mjy(O) is also presented. Finally, 
in Section 5, we analyze the resultant 5-L/electroweak hierarchy and show 
that it is of order 10 to 10^. The complete spectrum of squark and slepton 
masses, written in terms of the B-L boson mass, is then evaluated. Our anal- 
ysis depends on a numerical solution for the Higgs mass parameter, m/^(t)^. 
This is discussed and presented in Appendix A. In Appendix B, we numer- 
ically calculate the relationship between m^{Q) and m//(0) used in the text. 
Finally, in Appendix C we verify that our B-L / electroweak breaking vacuum 
satisfies the standard constraint and minimization equations presented, for 
example, in [601 ESI [M] . 
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2 The = 1 Supersymmetric Theory 

We will consider an = 1 supersymmetric theory with gauge group 

G = SUi3)c X SUi2)L X f/(l)y X f/(l)B_L (1) 

and the associated vector superfields. The gauge parameters are denoted 
by gs, g2, Qy and Qb-l respectively. The matter spectrum consists of three 
families of quark and lepton chiral superfields, each family with a right-handed 
neutrino. They transform under the gauge group in the standard manner as 

Q, = (3,2,1/3,1/3), n, = (3,1,-4/3,-1/3), = (3, 1, 2/3, -1/3) (2) 

for the left and right-handed quarks and 

L, = (1,2,-1,-1), z/, = (l,l,0,l), ei = (l,l,2,l) (3) 

for the left and right-handed leptons, where 2 = 1,2, 3. In addition, the spec- 
trum has one pair of Higgs-Higgs conjugate chiral superfields transforming 
as 

if =(1,2, 1,0), ^=(1,2,-1,0). (4) 
When necessary, the left-handed SU(2)l doublets will be written as 

Q, = iUi,D,), Li = {Ni,Ei), H={H+,H'), H={H\H-). (5) 

There are no other fields in the spectrum. 

The supersymmetric potential energy is given by the usual sum over the 
modulus squared of the F and D-terms. In principle, the F-terms are de- 
termined from the most general superpotential invariant under the gauge 
group, 

3 

W = fiHH + ^ [Xu^ijQiHuj + \d,ijQiHdj + X^^yLiHuj + Xe^ijLiHcj) (6) 

Note that an innocuous mixing term of the form LiH, as well as the dangerous 
lepton and baryon number violating interactions 

LiLjCk, LiQjdk, Uidjdk (7) 
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which generically would lead, for example, to rapid nucleon decay, are disal- 
lowed by the U{1)b~l gauge symmetry. To simplify the upcoming calcula- 
tions, we will assume that we are in a mass-diagonal basis where 

^u,ij = ^d,ij = = ^e,ij = 0, i ^ j- (8) 

Note that once these off-diagonal couplings vanish just below the compacti- 
fication scale, they will do so at all lower energy-momenta. We will denote 
the diagonal Yukawa couplings by = Xi, i = 1, 2, 3. 

Next, observe that a constant, field-independent /i parameter cannot arise 
in a supersymmetric string vacuum since the Higgs fields are zero modes. 
However, the HH bilinear can have higher-dimensional couplings to moduli 
through both holomorphic and non-holomorphic interactions in the superpo- 
tential and Kahler potential respectively. When moduli acquire VEVs due 
to non-perturbative effects, these can induce non-vanishing supersymmetric 
contributions to /i. A non-zero /x can also be generated by gaugino condensa- 
tion in the hidden sector. Why this induced /x-term should be small enough 
to be consistent with electroweak symmetry breaking is a difficult, model 
dependent problem. In this paper, we will not discuss this "/i-problem" , but 
simply assume that the /i parameter is at, or below, the electroweak scale. In 
fact, so as to emphasize the B -L /electroweak hierarchy and simplify the cal- 
culation, we will take fi, while non-zero, to be substantially smaller than the 
electroweak scale, making its effect sub-dominant. This can be implemented 
consistently throughout the entire scaling regime. The exact meaning of 
"sub- dominant" is quantified in Appendix C, where we also present the up- 
per bound on fi and, hence, the Higgsino mass in our approximation scheme. 

The SU{3)c and SU{2)l D-terms are of the standard form. We present 
the U{1)y and U{1)b-l -D-terms, 

DY=iY + gY(i>\{Y/2)^^ct>B (9) 

and 

Db~l = is^L + 9b~l<P\ (Yb~l)ab 0b (10) 

where the index A runs over all scalar fields (pA-, to set the notation for the 
hypercharge and B-L charge generators and to remind the reader that each of 
these D-terms potentially has a Fayet-Iliopoulos additive constant. However, 
as with the /x parameter, constant field-independent FI terms cannot occur 
in string vacua since the low energy fields are zero modes. Field-dependent 
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FI terms can occur in some contexts, see for example [OS]- However, since 
both the hypercharge and B-L gauge symmetries are anomaly free, such field- 
dependent FI terms are not generated in the supersymmetric effective theory. 
We include them in fl9|l. f|T0|l since they can, in principle, arise at a lower scale 
from radiative corrections once supersymmetry is softly broken [57] . Be that 
as it may, if calculations are done in the D- eliminated formalism, which we 
use in this paper, these FI parameters can be consistently absorbed into the 
definition of the soft scalar masses and their beta functions. Hence, we will 
no longer consider them. 

In addition to the supersymmetric potential, the Lagrangian density also 
contains explicit "soft" supersymmetry violating terms. These arise from the 
spontaneous breaking of supersymmetry in a hidden sector that has been inte- 
grated out of the theory. This breaking can occur in either F-terms, D-terms 
or both in the hidden sector. In this paper, for simplicity, we will restrict 
our discussion to soft supersymmetry breaking terms arising exclusively from 
F-terms. The form of these terms is well-known and, in the present context, 
given by HSl Sni EDI [Sa ED 

Vso{t = V2s + V3s + V2f, (11) 

where V2s are scalar mass terms 



1 i^l^Qil^ 
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(12) 

^h\H\'^ + ml\H\^) - {BHH + he), 
Vss are scalar cubic couplings 

3 

Vg. = Y^{A^^QiHui + Ad,QiHbi + A^^L.HDi + A.^L^Ha + he) (13) 

i=l 

and V2/ contains the gaugino mass terms 

V2f = ^MgAgAg + ^MsAaAa + ^MyAyAy + ^Mb-l\b-l\b-l + he. (14) 

As above, to simplify the calculation we assume the parameters in f|T2|) and 
f|T3|) are flavor-diagonal. This is consistent since once the off-diagonal pa- 
rameters vanish just below the compactiflcation scale, they will do so at all 
lower energy-momenta. 
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3 The Renormalization Group and B-L 



In this section, we discuss the spontaneous breakdown of the gauged B-L 
symmetry. The parameters in our theory all scale with energy-momentum, 
each obeying the associated renormalization group equation (RGE). In this 
section, we will solve those equations required in the analysis of B-L breaking 
to the one-loop level. 



Gauge Parameters: 

We begin by considering the RG running of the gauge coupling parame- 
ters. Since our low energy theory arises from an 50(10) compactification of 
heterotic string theory broken to SU{3)c x SU(2)l x U{1)y x U{1)b-l by 
Wilson lines, it is conventional to redefine the hypercharge and B-L gauge 
parameters as 

9i = yi^^' 9a = ^9b~l ■ (15) 

With these redefinitions, and defining t = In(^), the four running gauge 
parameters ga{t), a = 1, . . . , 4 all unify to a value g{0) at a scale M„. Precision 
measurements set [531 ESI EZ] 

g{0) ~ .726, M„ ~ 3 X lO^^GeV . (16) 

Note, g{0) is simply obtained from [23] and the equations there in. For 
specificity, we will use these values in our analysis, ignoring as sub-dominant 
the defocussing effects of possible string thresholds and non-universal soft 
breaking parameters. The RGEs are given by 

^ = J_f3^^ a = l,...,4 (17) 

where 

33 

(3a = bagl 6 =(-,1,-3, 12). (18) 
These can be integrated directly to yield 

9a{tr= a = l,...,4. (19) 
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In this section, we are interested in scaling all parameters from the unifi- 
cation mass Mu ~ 3 x lO^^GeV to the B-L breaking scale Mb-l — lO^Gel^; 
that is, 

tB^L ^ -28.7 <t<Q = tu. (20) 

In this range, we find from (|T8l) . ( |T9i) that the infrared free parameters gi^ g2 
and (^4 decrease for small energy-momentum as 

.441 < ^ < 1, .839 < ^ < 1, .303 < ^ < 1 (21) 

whereas the asymptotically free coupling g^ grows as 

2-36 < 1^ < 1 • (22) 

Gaugino Masses: 

Denoting My = Mi and Mb-l = M^, the RGEs for the gaugino masses in 
f lMj) are [55] 

^._L,„,2M., a.l,...,4 (23) 
where the ba coefficients are given in ffTSjl . This is immediately solved to give 

^^^^ = 7^^^' a = l,...A- (24) 

A priori, there is no constraint on the initial values Ma{0). In the scaling 
range ([201), it follows from ([MD that 

MUt) AUt) MAt) , , 

.441 < — if4 < 1, .839 < —44 < 1, -303 < — 44 < 1 25 

- Mi(0) - ' - M2(0) - ' - M4(0) - ^ ^ 

and 

2.36 < ^ < 1 . (26) 

- M3(0) - ^ ^ 

As will be seen shortly, the quantity we will be most interested in is 
ga{tf\Ma{t)\^ . Using ([IHD and ([24D, this combination runs under the RG as 
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Note that even if one assumes that the gaugino masses are "unified" at t = 0, 
making any ratio g°|o):a[j^°(Qjp unity, it is clear that the gluino mass contribu- 
tions will quickly grow to dominate. For example, at the electroweak scale the 
ratio of the gluino to the SU{2)l gaugino terms is 25.6. In this paper, so as 
to simplify the calculation and allow for a quasi-analytic solution, we will not 
assume unified gaugino masses, instead taking |Mi(0)p, |M2(0)p, |M4(0)p <S 
|M3(0)p. It then follows from ([27]) that 

gl\M^\\gl\M2\\gl\M,\^ <^gl\M:,\'' (28) 

over the entire scaling regime. Recall that "non- unified" gaugino masses 
easily occur in string vacua, while unification requires additional "minimal" 
criteria [60l[6l]. These are not generically satisfied in our MSSM theory. 



Slepton Masses: 

The RGEs for the slepton mass parameters m^,., mg. and my. are given by 



IStt^^ = 2{ml^+ml + ml;)\K^\^ + 2{ml^+m\ + ml)\K^\^ 
+2|A,J2 + 2|AeJ' - 2gl\M^\^ - Qgl\M2\^ - ^gl\M^? 
-\glS-\glS\ (29) 



~9l\MA' - Iglmi^ + IglS + IglS' , (30) 



■Iglmi^ + lglS' (31) 
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where 



o 2 2 , V^/ 2 n 2 , 2 2 , 



i=l 



Tri—m^ 



(32) 



5' 



f^l — iTT''^ — 2m r + rnt + 

Hi iJi c-i f 



i=l 



Tr{YB-Lm^ 



(33) 



A full numerical solution of these equations will be presented elsewhere. 
Here, we give an approximate solution based on the following observations. 
First, note that the initial conditions for the A-coefficients in equation f|T3!) 
are, ignoring phenomenological constraints for the time-being, completely 
arbitrary. However, it is conventional [M] to let 

An, = XuA.. = A., = XuA.. ^e, = Ae,ie, (34) 

where the dimensionful A-parameters satisfy 

Au, ~ 0(m„J, ^ ~ 0(mdj, A,^ ~ 0(m,J, i^, ~ 0{m,J . (35) 

This is not a requirement in the "non-minimal" string vacua that we are 
discussing. Be that as it may, for simplicity of presentation we will assume 
(!34|) and (!35|) for the remainder of this paper. Having done this, it follows 
that every term on the right hand side of equations (!29|) . (!30|) and (!3T|) . with 
the exception of the terms involving the gaugino masses, has the form of 
either |Apm^ or g'^im?. Our second observation is that the Yukawa couplings 
appearing in fl29l) . fl30|) and fl3Tl) satisfy 

|A^J < |A^J < |A^3| ~ 10"^ < ^a, |AeJ < lAeJ < lAeal ~ 10"^ < (36) 

throughout the scaling range fl20l) for a = 1, . . . , 4. Using flM|) . fl35|) and fl36|) . 
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it follows that one can approximate the slepton mass RGEs as 
levr^^ ^ -2glm' -Qgl\M,\' - lgl\M,\' 

-\glS-\glS' ^ (37) 
ISvr^^ ^ -^{9l\MA'-\glm^ + \glS + \glS' , (38) 

16vr^^ - -^^?I|M4r + ^^74^5'. (39) 

Third, recall that the initial gaugino masses M(i(0), a = 1, . . . , 4 are chosen 
so that (128|) is satisfied, but are otherwise arbitrary. Henceforth, we further 
restrict them so that 

gl\M,\', gl\M2\\ gl\M^\' « gjS' (40) 

over the entire scaling range ( l20l) . Fourth, we make a specific choice for the 
scalar masses at the unification scale M„. These are taken to be 

mniOf = m5(0)^ mg^Of = m^^{Of = m,,(0)2 (41) 

and 

mi,(0)^ = m,^(0)V^..(0)' (42) 
for all i,j, k = 1,2, 3. Note that the sneutrino masses are different than those 
of the remaining sleptons. This asymmetry is one ingredient in breaking 
U{1)b-l o,t an appropriate scale. Other than that, this choice is taken so 
as to simplify the RGEs as much as possible and to allow a quasi-analytic 
solution. We point out that soft scalar masses need not be "universal" in 
string theories, since they are not generically "minimal" . We emphasize that 
a 5-L/electroweak hierarchy is possible for a much wider range of initial 
parameters. 

Finally, let us consider the glS and gIS' terms. Note from (l32l) and (133!) 
that if one chooses the initial scalar masses to satisfy (14T|) and (l42l) then 

5(0) = 0, 5'(0) = Y.{~m,X^f + m,(0)^) ^ . (43) 

i=l 

Additionally, we will take 5' > with the scale set by the initial sneutrino 
masses. It then follows from ( H3l) that over the scaling range (l20l) 

glS « glS' . (44) 
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In fact, we can go one step further. Decompose S' in fl33p as 

S' = S'o + S'l, 

where 



3 



i=l 1=1 

Then we see from fl41l) and fl42l) that 



S'oiO) = 0, S\iO) = J2i-rnLX0r + m..m ^ 

i=l 

and, hence, over the scahng range fl2U]) 

qIS'q < 6-4 5' 1 . 

We conclude from (140!) . (144|) and (HHj) that a further approximation 
slepton mass RGEs is given by 



From (06]), (gnD and dSI]) one finds a RGB for S'l given by 



dS'i 9 2 c/ 



(it 327r2 
Using f[T^ . this is easily solved to give 

It follows from (fT6!) and (fT8|l that in the scaling range (l20l) 
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It is now straightforward to solve fl49p . fl50p and 0511) for ml, and m^, 
respectively. From (IT^ . (1501) and (1H5]) one finds 

= m,.(0)^ + 1(1 - (1 - ^^)-^/^^^)5',(0) (55) 

and 

mjt)^ = m,(0)^-l(l-(l-^^r/^^^)5\(0) (56) 

m..(t)^ = m.,(0)^-l(l-(l-^^)-^/^^^)5',(0). (57) 

These equations are almost identical with the exception of the sign of the 
second term. This is positive for m^^, whereas it is negative in the expressions 
for 777,2. _ This sign has important physical consequences. Using (fT6!) 

and flTHl) . the mass parameter parameter for m^^ increases from its initial 
value at t = to 

mLAtB-L? = trlM'' + (3.35 x 10-^)S\{0) (58) 

at ts-L- Thus it always remains positive. On the other hand, the mass 
parameters for ml. and ml. start at their initial values at t = but decrease 
to 

m,XtB-Ly = m,XOf - (3.35 x 10-^)S\{0) , (59) 
m^Xts^L? = m,^(0)2-(3.35x 10-2)5'i(0) (60) 

at ts^L- If the initial masses are sufficiently small and iS'i(O) is sufficiently 
large, then one or more of these squared masses can become negative signaling 
possible symmetry breaking. 

Spontaneous B-L Breaking: 

Given the RGE solutions in the previous subsections, one can now study 
the scalar field potential at any scale in the range (120|) and, in particular, 
search for local minima. Here, we limit the discussion to the slepton fields. 
Higgs fields and squarks will be discussed in the next subsection, where the 
consistency of this two-step procedure will be demonstrated. 
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Let us begin by considering the quadratic mass terms near the origin of 
field space. The relevant part of the scalar potential is 

V = V^, + ^Dl_^ (61) 

where and Db^l are given in (fT2|) and (fTOj) respectively. Expanding this 
using the B-L quantum numbers listed in (jl]) , the quadratic terms are given 
at any scale t by 

3 
i=l 

The reader should recall that the FI terms have been absorbed into the soft 
mass parameters. All (3 functions in the D-eliminated formalism are written 
in terms of these redefined masses. As the theory is scaled from t = 
toward t_B_L, the m|., ml. and ml. parameters scale as in (!58l) . (!59l) and (!60l) 
respectively. 

The first requirement for spontaneous B-L breaking is that at least one 
of the slepton squared masses becomes negative at t^-L. Clearly, this can- 
not happen for mi^.{tB-LY ■, which is always positive. However, if the ini- 
tial squared masses are sufficiently small and iS'i(O) sufficiently large, both 
me.(tB-L)^ and my-itB-iY can become negative. Since the fields are elec- 
trically charged, we do not want them to get a vacuum expectation value 
and, hence, we want me-(tB-L)^ to be positive. On the other hand, the Vi 
fields are neutral in all quantum numbers except B-L. Hence, if they get 
a nonzero VEV this will spontaneously break B-L at Ib-l-, but leave the 
SU{?>)c X SU{2)l X U{1)y gauge symmetry unbroken. We now show that 
for a wide range of initial parameters this is indeed possible. For simplicity, 
let us choose the initial right-handed slepton masses to be 

m,,(0) = m,,(0) = Cm,(0), 771,3(0) = m,(0) (63) 

and 

m,M = rrieM = m,,{0) = Am,(0) (64) 

for some dimensionless constants C and A to be determined. Using P2|) . 
fH3|l and fl63|l . flMl) we see that S'i{0) is parameterized by 

S'liO) = {1 + 2C^ -3A^)m,{0y . (65) 
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Let us first consider 1/3. It follows from fl60p . 0631) and that 

mu^{tB-Lf = (1 - (3.35 X 10-2)(1 + _ ^A^))m^{Qif . (66) 
For specificity, we will henceforth take 

(3.35X 10-2)(l + 2^2-3^2) = 5 . (67) 
This choice yields the simple result that 

m,^{tB^Lf = -^mM\ (68) 

suggesting a non-zero VEV in the 1/3 direction. Now consider z/i, V2- In this 
case, we see from fl60l) . fl63l) . fl65|) and fl67j) that 

m,,,{tB-Lf = (C' - 5)m,(0)2 . (69) 

Now let us examine the 61^2,3 fields. Using f lS^ . flMl) . fl^ and flB7|) one finds 

me,(tB_L)' = (^'-5)m,(0)2 (70) 

for i = 1,2,3. Since all me-{tB-L)^ must be positive, the coefficient A must 
satisfy — 5 > 0. Again, for specificity we will choose 

A = V6 , (71) 

which yields the simple result that 

meXtB^L? = m,{Of (72) 

for i = 1, 2, 3. Putting A = \/6 into expression (1671) gives 

C = 9.12 . (73) 

Hence, we see from f l69l) that both my^^{tB-L)'^ are positive and given by 

m,,,{tB^Lf = l?>.2m,{Qf . (74) 

We conclude from ( 172|) and ( 174|) that, near the origin of field space, there are 
positive quadratic mass terms in the 61^2,3 and 1/1^2 field directions. Finally, 
let us consider the L^, i = 1, 2, 3 masses.' Using i^, ([58]), dSSD, ([67D 

and (17T|) . we find that 

mL,,,^,{tB-L)'' = II m,{0)\ (75) 
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that is, a positive quadratic mass term in each if the Lj field directions at 
the origin. We note for future reference that using (17T]) and (175]) . equation 
( 165|) becomes 

S'i{0) = 149 m^(0)2 . (76) 

Given these mass terms, as well as in (12T1) . one can minimize the com- 
plete potential to determine the vacuum at ts-L- The part of this potential 
relevant to finding a local minimum in the slepton fields is V in (I^Tj) . We note 
that, in principle, the |Ai/3p|L3p|i^3p term in the |(9hW^P contribution to the 
potential, where W is given in ([6]), could play a role in selecting the vacuum. 
However, since X^.^^ <^ g^, this term can be safely ignored. Expanding out the 
contributions to V and collecting terms, the important interactions are 



T/ 2 I |2 I 4,94 I |4 I / 2 I 2 1 |2m |2 

y = mjusl +^|z/3| +(m^,^, + -^4^31 )ki,2| 



+ E(K + l9>3\')\e.\' + - IglHW) + ■ ■■ (77) 



i=l 



Recalling from (168!) that m^^itB^LY = — 4mj,(0)^, the first two terms in (1771) 
can be minimized by the non-zero VEV 

= . (78) 



Is this point a local minimum of the complete slepton potential? To deter- 
mine this, consider potential (1771) near this VEV. Clearly, the first derivatives 
of all fields vanish at this point. Now determine the masses in all slepton 
field directions evaluated at (178|) . It follows from (177|) that 

V = + {mlJM' + + iO\L:\') + • • • , (79) 

where 



1=1 



3 



«) = ml + lgl{u,r, {ml)=ml-\gl{v,f . (80) 
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Using dM]), (IZ2D, dUD, dZS]) and dZHD we find that at tB^i 



= 8m,(0)2, {mlj = 82.2 mM', 
{ml) = 5m,(0)2, (miJ = 7m,(0)^ (81) 

Since all of these masses are positive, we conclude that the vacuum specified 
by 

M = 0, (us) = ^ ^ (82) 



B-L 



4 

and 

(e,) = (Li) = 2=1,2,3 (83) 

is indeed a local minimum of the slepton potential V . This vacuum sponta- 
neously breaks the gauged B — L symmetry while preserving the remaining 
SU{?>)c X SU{2)l X U{1)y gauge group. The massless Goldstone boson is 
"eaten" by the B-L vector boson giving it a mass 

Ma^_^ = V2 gB^LitB^L) {i^z) . (84) 

Using (fTSjl and (182|) . this becomes 

MA^_,=2V2m,{Q) . (85) 

We now have to include the Higgs fields and squarks, analyze their masses 
at tfi-L around vacuum (182!) . (!83|) and show that this is a local minimum in 
the complete scalar field space. To do this, one must discuss the Higgs and 
squark RGEs, to which we now turn. 



4 The RGEs for Higgs Fields and Squarks 

In this section, we present and analyze the RGEs for the up and down Higgs 
fields, the left and right squarks, as well as the associated fi and B parameters. 
Using these results, we compute the Higgs and squark masses at t^-L and 
show that they are all positive at vacuum fl82l) . fl83|) . as desired. 
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"Up" Higgs and Squark Masses: 

The RGEs for the "up" Higgs and squark mass parameters uih, tuq^ and 
m„- are given by [55] 



lQ^2d^ = E(6K,+< + 0|A..r + 2K^+m^ + m2j|A,j2 

i=l 

+6\A^f + 2\A,f) - ^gl\M^\'-Qgl\M2\' (86) 

+2|A„J2 + 2\A,f - ^^^glm\' - 6glm' - y^^a'lMap 
-l9lm\' + ^glS + \glS' , (87) 
levr^^ = ^ml^+ml + ml)\X^f + 4\A,f 

-fgim' -Igim' m 

—gls - ^gls' 

where S and S' are given in ( l32l) and ( l33l) respectively. 

A full numerical solution of these equations will be presented elsewhere. 
Here, we give an approximate solution based on the following observations. 
First, as discussed in the previous section we will assume that the A-coefiicients 
satisfy (IMI) and (135|) . Having done this, it follows that every term on the 
right hand side of equations (!86l) . (IHTl) and (IHHI) . with the exception of the 
terms involving the gaugino masses, has the form of either |Apm^ or g'^rn^. 
Our second observation is that the Yukawa couplings appearing in ( 1861) . (IHj 
and dHH]) satisfy 





< \Xu2\ 










< -^^2 


< 


~ 5 X 10~2 < ^a, 


(89) 



and 

|A„J < |A„,| ~ 10-2 < |A„J ~ 1 (90) 
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throughout the scahng range (120]) for a = 1, ... ,4. We see from (Mli) that 
lA^gl is o/ (9(1) and, hence, terms containing it cannot be dropped from the 
RGEs. As in the previous section, we will continue to assume that the 
initial gaugino masses M(j(0),a = 1,2,4 are chosen sufficiently small that 
the inequalities in PUI) remain satisfied. However, we will not make a similar 
assumption about M^(fi), allowing it, for the time being, to be arbitrarily large 
and not sub-leading to S'l . Since the initial scalar masses satisfy fHTl) and 
(142|) . the inequalities (jHj) and (HHj) are still satisfied. Finally, we will assume 
that terms with gfS are small compared to any term proportional to lA^jgp. 
Using all these inputs and assumptions, it follows that one can approximate 
the squark mass RGEs as 

levr^^ ^ 2{ml,^+ml + ml)\X^.f + 2\K.f\A^f 



dt 



+ (91) 



-f9l\Ms\'-\9lS\ (92) 



and 



16-^^ - -f9lm\' + \9lS\, (93) 

16-^^ - -Y9lm'-\9lS\. (94) 

Let us begin by analyzing equations fl9T]) and fl92|) . Note that fl9T]) and 
fl02]) can be written as 

16,=_J?i ^ 16-=^ - jsm? + J*V. . (95) 

^ (96, 

respectively. These are easily solved to give 



1 9 

f^n^ — —rriTT ;r 

m^, ~ -m%, 



1 + 1^ glS\ + Cq, , (97) 
£ 9l\Ms\' 1 9lS\ + (98) 
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where 



* 2| - ' '-^ 8 , 1 



are evaluated using f|T9|) . fl2^ and fl53l) . The integration constants are 

CQ3=mQ3(0)2-lmH(0)2, C„3 = mQ3(0)2-^mH(0)2 (101) 



where we have used our assumption fHTj) that the initial squark masses are 
degenerate. Note that since 63 = —3 and t < 0, expression fl99|) is always non- 
negative. Similarly, since 64 = 12, it follows that fllOOl) is also non-negative. 

Now consider the m'jj equation flS^ . If we insert fl^ and fl^ into fISB]) . 
we get an expression for without the the squark mass squared terms. 
It is important to note that in doing so, we would gain a term in the beta 
function that depends on the gluino mass squared. If we then apply (HOjl and 
the discussion above (PT]) (where we took the gluino mass squared term to 
be of the same order as S'l in our approximation), we get 

A|A„3P(2m^-^ r,2|^3|2 + ^2 + |i |2) (102) 
at ovr^ dvr^ Jo 

with rri^ defined by 

ml = Cq, + Cu, = 2mQ3(0)2 - mH^Q? . (103) 
Henceforth, to simplify our discussion, we assume 

rn,M' = . (104) 

which sets 

ml = . (105) 
Furthermore, we will choose the initial value ^1,13(0) so that 

|i„3p « 2ml - ^ fglmi' (106) 
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over the entire scaling range. As a final simplification, we note that A^g scales 
by only a few percent over the scaling range. Hence, we approximate it as a 
constant with its phenomenological value of 

A„3 - 1 ■ (107) 
Using (fT05D . (fT06D and (fTOTI) . RGE equation (fT02|) simplifies to 

^-ip"i,-^r«. (108) 



In this paper, we will solve equation (11081) subject to the constraint that 
rrijj, which is positive at t = 0, remain positive over the entire scaling range. 
Then, RGE (11081) is equivalent to an integral equation for given by 



_ 3 rO(^,,-^fo silMal^ 

ml c^mniOYe ^ . (109) 



Note that (HM^ imphes 

Cqs = IrnniOf, C^, = -^"^^(0)^ . (110) 
It follows that the m^^ and m^^ equations in fPTI) and (jHSD become 

respectively. Finally, let us consider RGBs (p3l) and These are easily 

solved to give 



- ^3^3^ + ^ giS\+'^^^ (113) 



where the (yf||M3p and g\S'i integrals are given in fl99l) and (llOOp respec- 
tively and we have used assumption (l^Tl) that the initial squark masses are 
degenerate. 
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"Down" Higgs and Squark Masses: 

The RGEs for the "down" Higgs and squark mass parameters mg and md. 
are given by [55] 



IGtt'^ = ^(6(m^^+m^ + mi)|A,j2 + 2K^+m| + m2j|Ae,r 

i=l 

+Q\A,f + 2\A,f) - -gl\M,\'' - GgllM^l' (115) 

-^gfmi' - f ^l|M3r - lalm' (ne) 

where S and S' are given in (l32l) and (133|) respectively. First consider the 
rud^ equation. Using the assumptions hsted in (IMI) . (!35|) . (HOj) . (l48l) and (!89|l . 
equation ( I116P can be approximated by 

ISvr^M ^ --^73^31' - -9lS\ . (117) 
This can be immediately integrated to give 

where the first and second terms on the right hand side are evaluated in 
(l99l) and fllOOl) respectively, and we have used assumptions (jHj) and fll04p . 
Note that these are both non-negative and, hence, m'^,,i = 1,2,3 are all 
non- negative throughout the entire scaling range (120|) . Now consider the rriff 
equation flll5p . All of the terms on the right hand side of this expression are 
small compared to (7IIM3P and gfS'i. In addition, these terms are multiplied 
by which is small. Hence, to the order we are working 

—^-0 (119) 
dt ^ ' 
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and 

where we have used (HTl). 



|~mH(0)^ (120) 



m 



The fi Parameter: 

The fi parameter enters the potential for the Higgs supermultiplets H, H 
through the superpotential in ([6]). The RGE is given by [55] 

This is easily integrated to 

= ^(o)e^Ti^^"(^-i(=^l^-»l'+=^l^'*.l'+l^^'l'+l^=«l')-t^??-3<;i) _ (122) 

As discussed earlier, a constant field-independent /i(0) parameter cannot arise 
in a supersymmetric string vacuum since the Higgs fields are zero modes. 
However, a yU-term can arise from higher- dimensional couplings to moduli. 
Although typically chosen to be of electroweak order, in this paper, to em- 
phasize the 5-L/electroweak hierarchy and simplify the calculation, we will 
choose yu(0) to be substantially smaller, making its effects sub-dominant. It 
is clear from fll22p that fi runs slowly over the scaling range and, thus, once 
chosen to be sub-dominant at the unification scale it remains so throughout 
the scaling range. The exact meaning of the term "sub- dominant" will be 
given in Appendix C. 

The B Parameter: 

The B parameter enters the potential for the Higgs scalars H, H through 
the soft supersymmetry breaking term V2s in (fT2|) . The RGE for B is given 
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by m 



16-^^ = B(j2mx^f + ^\>^^f + \Kf + \>^ef)-U-^9i 

\i=l 

+ Jy.^6A^^^:^ + A,X.+^KK+^A,XJ (123) 
V i=i 

+ 6glM2 + ^glMi 

Recalling the discussion in the previous subsection, we take the term propor- 
tional to /i to be sub-leading. Equation (11231) is then solved by 

B = S(0)e-Ti^^°(^-i('l^-'l'+='l^''-l'+l^^'l'+l^'='l'^"^^'-'^') . (124) 

Here, unlike the constant, field-independent /i(0) parameter, the dimension 
two 5(0) parameter arises from supersymmetry breaking and need not van- 
ish. Using (1361) . (I89p and (1901) . the RGE for B can be approximated by 

B ~ B(0)e"T^^°l^"3l' . (125) 



This parameter will make its appearance in the analysis of the Higgs potential 
below. The assumption that the term proportional to /i is sub-leading is 
easily checked using (I106P and the relative sizes of B and fi presented in 
Appendix C. 



Higgs and Squark Masses at the B-L Breaking Vacuum: 

At the end of the previous section we showed that, when restricted to slepton 
scalars only, the potential energy has a local minimum given in (182|) and (!83|l . 
Clearly, however, to understand the stability of this minimum it is essential 
to extend this analysis to the entire field space; that is, to include all Higgs 
and squark scalars and as well as the sleptons. Given the RGE solutions in 
the previous subsections, one can now study the full scalar potential for any 
scale in the range fl20l) . 

Let us begin by considering the quadratic mass terms, not at the origin 
of field space but, rather, at the slepton minimum given in fl82|) and fl83|) . 
The relevant part of the scalar potential is still flMl) . However, this is now 
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evaluated for all scalar fields near the slepton VEVs. Note that, in principle, 
the |Ai/3p|ifp|z/3p term in the |c}2,3iyp contribution to the potential, where 
W is given in ([6]), could play a role in selecting the vacuum. However, since 
■^U3 <^ 94:, this term can safely be ignored. Similarly, under the assumptions 
flM|) and fl35l) . one can neglect the A^^L^Hu^ + he term in V^s- Expanding 
out V using the B-L quantum numbers listed in (jlj), the quadratic terms at 
any scale t are 

V = --- + V^2 +V^2 +V^2 ^+..., (126) 

slepton Jrliggs squarks 

where IC,2 is given in (1791) . (IHOl) . The Higgs contribution to the quadratic 

slepton 

potential is 

V^2^,^^^ = m]j\H\^ + m\\H\^ - B{HH + he) (127) 

where mj^, and B are given in fll09p . (11201) and (I125P respectively. Note 
from fll09p . (11201) that for t <C the quantity |m|^ — m|^| becomes non-zero 
and large. Henceforth, we will assume that in this range of t the coefficient 
B is such that 



B 



1 1 



<1. (128) 



This is easily arranged by adjusting B{0). For t ^ 0, the Higgs mass matrix 
(11271) can be diagonalized to 



Kn2 =mUH'\' + mji,\H'\\ (129) 



where 



m\, c:i m\ — m%r ( — ^ ^) , m L ~ m f> — m?r ( — ^ ^) (130) 



and 



H' c^H+[ , ^ A h*, H'-( „ ^ , ]H* + H . (131) 

H ~ H J \ ~ H - 



Here, H*^ = g^^^e^nH*^ and H*"^ = g^^ecoH*'^ have the same SU{2)l x 
U{1)y transformations as H and H respectively. In component fields 

H* = {H-\ -ffO*), H* = (i/°*, -//+*) . (132) 
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It follows from ffTU^ . ffT^ . ffT^ and ffMID that anywhere in the range 
t -C 0, and specifically at ts-L, 



mjj, :^ mjj = mniO) > ■ (133) 

Importantly, however, we see from (11091) . (I120p that as t becomes more neg- 
ative can approach, become equal to and finally become smaller than 
m^^ {B / {m^g — m\)Y ■ This corresponds to rn^j, being positive, zero and 
negative respectively. The vanishing of m|^, signals the onset of electroweak 
symmetry breaking. Clearly, to evaluate m?^, at Ib-l we must solve equa- 
tion (llOSp or, equivalently, (11091) for the running of Ta\. The solution of 
this equation will depend on two arbitrary parameters, m//(0)^ and, using 
^^^1 M3(0). It follows from the exponential form of (I109p and the fact that 
integral (!99|) is non-negative that increasing |M3(0)| decreases muitY for 
any fixed value of t, and vice versa. Hence, specifying |M3(0)| is equivalent 
to specifying the value of niHitY at some fixed t. For reasons discussed in 
in detail below, the physics is most transparent if we choose this to be the 
electroweak scale corresponding to Iew — —33.3. Specifically, we take 

mnitEw? = (134) 



where 



T is defined by 



= • (135) 



^ ( , ) ~ 40 , (136) 

thus explicitly satisfying constraint (11281) . A is a parameter with the range 
< < 1. The upper bound follows from our requirement that m'jj be 
positive over the entire physical scaling range and, hence, positive at tEw- 
Note from ([TSOD, ([T33D, ([IMD, (^M and (^M) that 

-A^ 

It follows that electroweak symmetry will be broken at Iew only if A^ is 
strictly positive. Hence, the lower bound on this parameter. 

The role of (I134p in electroweak breaking will be thoroughly discussed in 
the next subsections. Its relationship to |M3(0)| and, hence, to the solution 
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of (11081) for rriHitY will be derived in Appendix A. Here, we will simply state 
the result that specifying this value at Iew is equivalent to choosing 

|M3(0)|2 = . 0352(1 -H^)m^,(0)^ (138) 

Using this expression for |M3(0)|, we can solve (11081) numerically for any 
fixed value of T'. The numerical results can then be fit to a smooth curve. 
We find that the numerical data is well-represented over the entire scaling 
regime by 

m^itr = (1 - (1 - ' (139) 

where b is a function of T' of the form 

b(r) = 19.9(1 - ^ '^^^ ) . (140) 

Note that at t = and t = tsw, ^'h is given by m//(0)^ and niHiOy /T''^ 
respectively, as it must be. One can now use expression (11391) to evaluate 
and, hence, m|^, at ts^L for any T and parameter A^. It suffices here to give 
the lower bound. To do this, first note that for = 1 the T dependence 
drops out of (I139p and (I140p . Furthermore, it follows from (11371) that this 
corresponds to the largest negative H'^ squared mass at tEw and, hence, a 
T-dependent lower bound on Tnjj, at ts-L- Taking A^ = 1, we find 



and, hence, 



mnitB^L)^ = .0565 mniOY (141) 
17.7. 



mH'itB-Lf = .0565(1 - ^) mniOf . (142) 

Since, from ([ISSD, ~ 40, this expression is always positive. We conclude 
that at tB-L, for any T and < A^ < 1, 

17 7 

mjj, > .0565(1 - — -) mnioy > . (143) 

Finally, note that evaluated at the slepton vacuum (152]) . (155]) the diagonalized 
quadratic Higgs potential remains 

=(m|,)|/7f + (m|,)|^f , (144) 
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with 

(m^,) = mjj,, {m%) = m% . (145) 
It then follows from fll43p and fll33p respectively that at Ib-l 

{m%) > .0565(1 - ^) mniOy, {ml,) = m^,(0)2 (146) 

and that they are both positive. 

The squark contribution to the quadratic potential is 



= + + , (147) 

4 = 1 

with 

(ml^) = + ^gl{u,)' (148) 

and 

/ 2 \ 2 1 2/ \2 

«) = rnl-\gl{u,)\ (149) 

Here mg . are given in f illip , (11 130 , . are given in f lll2p , (11140 and m^, are 
given in (IllSp . 

Using this result, one can now evaluate the squark masses in (I148P and 
(I149p . Let us begin by computing the second term on the right hand side of 
these equations. Using (fT6l) . (12T]) . (1761) and (fTSI) . we find that at t^-L 

\9>3)' = ^ m^iOY . (150) 
First consider ("ig^^)- Recall from (I113p that 

,2 2 /•* o,,.,2 , 1 2c/ , ^^(0)' 



rrin — 

VI, 2 



2 r<,,=iM3i=+ 1 r<,j5^+!^^ (151) 

For any value of t, the first term can be evaluated using fl^ . (11380 and the 
second with (^E>, fflOOD - We find that at te-L 

I ' ' 9lm\' = .143(1 - H;|)m^(0)^ (152) 
^JJ' alS'i = -1.66 m,(0)2 (153) 
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and, hence, 

2 57 

"^Qm - -643(1 - -Y^)niH{'^f - 1.66 m,{Qf . (154) 

Using this and ( I150p in 0148^ gives 

2 57 

(m^^ J = .643(1 - -^)mH{Qf - .336 mM' ■ (155) 

To further evaluate this squared mass and to determine whether or not it 
is positive, one must give an exphcit relationship between the parameters 
mj^(O)^ and mui^Y- This will be discussed in detail in Appendix B. Here, 
we will simply use the result. By demanding that all squark mass squares 
remain minimally positive at any scale Iew < t < 0, we find that 

2 25 

m^(0)2 = .864(1 - ^)mH(0)2 . (156) 
Using this to eliminate the m/f(0)^ parameter, expression fll55p becomes 

Km) = -408(1 -|5)m.(0)^ (157) 

Note that this expression is positive for all parameters ~ 40 and < 
< 1. Also, for sake of comparison to the sneutrino VEV, we will write all 
masses at ts-L in terms of m,^(0)^. The form of m^^ ^ ^ind m^^ for i = 1, 2, 3 
given in (11141) and (11181) respectively are similar to that of m^^ ^ and can be 
evaluated using (11521) . (11531) and (I156p . Using these results and (I150p . we 
find from (fT49|) that 

993 

«,) = K) = 1-08(1 - -^)mM ■ (158) 

These mass squares are positive for all ~ 40 and < < 1. 
Now consider {rnqj. Recall from (II lip that 

"L-^™?,- J./'9^|M3r + ^/ftV, + im„(0)= , (159) 

The second and third terms are given in (I152p and (I153P respectively. Unlike 
the masses evaluated above, however, ttIq^ also depends on the mass param- 
eter m'jj. This can be evaluated using (11390 . (11400 for any value of t. At 
ts-L "we find 

ml ^ .0565(1 - ^ + f^)mH{Or ■ (160) 
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Using this together with fll52p and (11531) . niQ,^ at t = ts-L is given by 

m^3 ^ -1.66m.(0)^ + 0.328(1 - - 1^)^^(0)2 , (161) 

which, using (11561) . becomes 

2 . 0.00185 0.545, ,^,2 

771^3 ^ -1.29(1 + + ^)m.(0)2 . (162) 

Similarly, m^^ given by (I112p is found to be 

2 , 0.00284 0.691, ,^,2 /..on 
m^3 ^ 1.68(1 - ^)^.(O)' (163) 

at t = tB~L- Using these results and (I150p . we find from (I148p . (ll49p that 

, 2 ^ . 0.0549 16.2, , ,2 

(m^3) ^ 0.0435(1 - - ^)m,(0)2, 

, 2 ^ / 0.0136 3.30, , ,2 , , 

(ml) ^ 0.353(1 - - ^)m,(0)^ (164) 

Note that they are both positive for all ~ 40 and < < 1. Putting 
everything together, the squark masses, expressed in terms of m,y(0)^, are 
given in (I157p . (11580 and (I164p and are all positive. 

From §2), ffT^.ffm. ffTFTD . ffT^ and ffTMD we conclude that at t = 
tB-L the vacuum specified by 

K2)=0, (z/3) = ^^^, (e,) = (L,) = (165) 

(Qi) = (ni) = {di) = (166) 

and 

{H') = {H') = (167) 

is a local minimum of the potential energy. It has positive mass squares 
in every field direction including H' and H' , signaling that, at t = Ib-l-, 
electroweak symmetry has not yet been spontaneously broken. 
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2 

Higgs 



Electroweak Symmetry Breaking: 

To explore the breaking of electroweak symmetry, one must add to 
in fll27p all other relevant interactions involving the Higgs fields. It follows 
from previous discussions that the relevant part of the Higgs potential is 

^ = ^'nW + ^^^ + ^E^W (168) 

a=l 

where, written in the mass diagonal fields H' and H' defined in ( 11311) . 



Dy = f~{\Hr-\Hr), (169) 

Dsu{2),a = ^{H'^cT.H' + H'^aaH') (170) 

and we have dropped terms of 0{T~^). The parameters in this potential 
should now be evaluated, not at ts-L, but, rather, at the electroweak scale 
Mew — lO^GeV which corresponds to a tsw — —33.3. The gauge couplings 
are easily evaluted there using (fT8|) and (fT9l) and found to be 

i^ = A05, ^%^ = .818. (171) 
^(0)^ ^(0)^ 

Also, it follows from ffT20D and ffWD that 

rriH'itEwf ^ rnniOf, mH'{tEwf = -^?ti/^(0)^ (172) 

where ~ 40 and < < 1. 

We can now proceed to minimize potential f ll68p . Writing 

H' = {H'+, if'°), H' = {H'\ H'-) , (173) 

potential (11681) becomes 

V = -^mH{0)\\Hy + \H'^\^) + mH{0)\\Hy + \H'-\^) 
1 ^ 







+ 


1(5 


8^5 






+ 




+ 


4\H' 
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This is easily minimized to give 



{{H")) = ^^ ^^^(0) ((i/'+)) = 0, (175) 

T ^Jlgi{tEwY + g2itEwY 

and 

((^'°)) = ((^'-)) = 0, (176) 

where the double bracket (( )) indicates the vacuum at tEw- K is straight- 
forward to compute the squared masses of the radial fluctuation 5H'^ and 
complex H'" fields at this vacuum. We find that 

V = {{m%,))\6Hy + {{m%,))\Hy + ((m|,.))|^'"|^ + . . . (177) 

where 

((,4,„)) = 4^« (178) 



and 



A 2 A 2 2 2 

(«,o)) = (l-:^)mH(0)^ {{m%^)) = il-^i^-\^))mHiOr . (179) 

Note that the Higgs mass squares are positive for all ~ 40 and < < 1. 
Evaluated at this minimum, the phase of if'" and H'~ , H'~* are Goldstone 
bosons which are "eaten" by the Higgs mechanism to give mass to the Z and 
vector bosons. For example, the Z mass is given by 

Although the mass eigenstate basis H', H' is the most natural for ana- 
lyzing this vacuum, it is of some interest to express it in terms of the original 
H and H' fields. Using p, and (dSn]), we find 

((iJ+)) = ((^-)) = (181) 

and, to leading order, that 

((^°)) = ^^"^^(0) m) = h{H^)). (182) 

T^Jlgi{tEw? + g2{tEw? ^ 
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Note that the condition {{H'^)) = in fll76p does not imply the vanishing of 
{{H^)). Rather, {{H^)) is non-zero and related to {(H^)) through the ratio 

|||^ = tan/3 = r + O(r-i) . (183) 

We have indicated the 0(T~^) contribution to emphasize that although 
tan/] = T to leading order, this relationship breaks down at higher order 
in T~^. We refer the reader to Appendix C for a more detailed discussion of 
this point. 

We conclude that electroweak symmetry is spontaneously broken at scale 
^Ew by the non- vanishing H'^ vacuum expectation value in (11751) . This 
vacuum has a non- vanishing value of tanj3 which, using the assumption for 
7"2 gjygj^ f|i3(j|) . satisfies 

tanf3 ~ 6.32 . (184) 

As far as the Higgs fields are concerned, the vacuum specified in (I175P and 
(11761) is a stable local minimum. As a check on our result, choose //^ of order 
0{T~'^) or smaller, that is, non-vanishing but sub-dominant in all equations. 
Then (11331) . (I134p . (11351) and (I136P satisfy the constraint equations, given, for 
example, in [61], for the Higgs potential to be bounded below and have a 
negative squared mass at the origin. Furthermore, to the order in we 
are working, (I180p and (11821) for the Higgs vacuum satisfy the minimization 
conditions in [51]. This is shown in detail in Appendix C. 



Slepton and Squark Masses at the EW Breaking Vac- 
uum: 

Clearly, to understand the complete stability of this minimum, it is essential 
to extend this analysis to the entire field space; that is, to include all slepton 
and squark scalars as well as the Higgs fields. The relevant part of the 
potential energy is 

V = V2s + \dI_^ + 1^2, + 1 ^ ^2^^^^^^ ^ (185) 

a=l 

where V2s, Db~l, Dy are given in f|T2|) . ffTOj) . Qj respectively and Dsu{2)La is 
the extension of (I170p to include all slepton and squark doublets. Note that, 
in principle, the superpotential can also contribute to the squark/slepton 
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potential energy. However, it follows from 0361) and fl89p . fl90p that such terms 
will be negligibly small compared to those in fll85p with the notable exception 
of the third family up-squarks. We will take this contribution into account 
when it arises below. Now expand V to quadratic order in the fields using 
the quantum numbers listed in ([2]), ([3]) and (jlj). We find that, evaluated at 
the slepton and Higgs VEVs in ([12]), ([HSD and ffT75]) . ffT75]) respectively, 

V = --- + V^2 +V^2 +V^2 ^ + . . . , (186) 

H iggs sleptons squarks 

where 14,2^^ ^ is given in ffTTTD . ffT78|) and ffTTOD . 

The slepton contribution to the quadratic potential is 

= (K3))i^^3p+(«,))i^i,2p+x:((K))i^^i' 

i=l 

H{rn%))m' + {{m%))\E,\') + . . . , (187) 

where 

(«)) = «), (KJ) = K,), 

(K)) = {<) + \^)9l{{H''))\ (188) 
(«)) = {^l)-\{l9l + 9l){{H''))\ 
(«)) = {0-\{\9l-9l){{H")? . 

The squared masses {ml,^), {ml^^), {ml,) and were defined in fISUl) and 

evaluated at ts-L in (181]) . Now, however, these values must be corrected 
by scaling down to tEw- Using (fT6|) and (fT8|) . the slepton masses defined 
in (!55|) - (!57|) can be evaluated at tEw- Note that the parameters A, C and, 
hence, iS'i(O) given in (17T]) . (1751) and (1751) remain the same. We find that 

ml^ = -4.38m,(0)2, ml^^ = 77.8m,(0)^ (189) 
ml = 0.625m^(0)^ m^^ = 11.4m^(0)^ 

Using these values, slepton potential (177|) has a local minimum at 

2m ((]) 

((z/3)) = (1.05) , ((^1,2)) = (190) 

j94itEw) 
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((e,)) = ((L,)) = ^ = 1,2,3 (191) 

We can now evaluate the squared slepton masses in (IHOl) at tEw- We find 
(m^J = 8.75m,(0)2, {mlj = 82.2m,(0)2, 
(ml) = 5m,(0)^ (ml) = 7 m,(0)2 . (192) 

Inserting these results into fllSSp and using fll75p gives 

{{ml.)) = 8.75m,(0)2, ((m^ )) = 82.2 m,(0)2, 



l^S' ' -f ' \\ 1^1,2/ 

(K)) = 5m.(0)^ + 2f^if^)^!l^, (193) 



{{m%)) = 7m,(0)2-- 



^2 ' 



{{mi)) = rmM' 



l9l + 9l ) 



Finally, using expression fll56p which relates m,^(0)^ to m/^(0)^ and fll7ip . 
these squared masses become 



{{ml;)) = 8.75m.(0)^ ((m^^ ^)) = 82.2 m,(0)^ 

{{ml)) = 5m,(0)2(l + A2^^:^)^5m,(0)2, (194) 

n 1 23 

{{m%)) = 7m.(0)2(l - A^^) ^ 7m.(0)^ 
{{m%)) = 7m.(0)^(l + A^^)^ 7m.(0)^ 

Note that the slepton squared masses are positive for all ~ 40 and < 
A2 < 1. 

Similarly, the squark contribution to the quadratic potential is 
Kn^Va.. = E((«))|f^^l' + («))lAP + («))|n.P + ((m^J)|d,p) (195) 



squark 

1=1 



with 



(«)) = {rnl)+[^^^ + \K\\s]{{H''))' 
{{ml)) = {ml)+(i^^){{H''))' (196) 
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and 



(«)) = {<)-{^^9l-\K?5^^{{H''))^ 

(«)) = K) + ^^7?((^'°))^ (197) 

Note that in the expressions for {{mfjj) and ((rra^^)) we have included the 
non-neghgible superpotential contribution. The squared masses {rriQ,), (m^.) 
and (m^.) were defined in ( 11481) . fll49p and evaluated at ts_L in (I164p . Now, 
however, these values must be corrected by scaling down to tEw- First, we 
must compute (74 at Iew- Using f|T6l) and f|T9l) we find 

'-^^ = .272 . (198) 

One can now evaluate the second term on the right hand side of (11481) and 
f lM using mm . The result is 

^glius)' = 1.46m,(0)2 . (199) 

The squark masses defined in (fTTTD. (fTT2|) . (fTT3|) .(fTT4D and ffTT8|) can be eval- 
uated at tEw using (I95]),(l96]),(l93l),(l9l,([in2l), ffTTSl) and ffT56l) . Adding these 
to ffT99ll . expressions ffT48ll . ffll9ll become 

{ml) = 0.132(1 - ^)m^,(0)^ 



4.87. 



^ J = 0.465(1 - ^)m^(O)^ 
7 72 

(m^3) = 0.374(1 - ^)m^,(0)^ (200) 

3 42 

{mlj = 1.041(1 -^)m^(0)^ 

3 49 

{ml) = 1.041(1 -^)m^(0)2. 
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Inserting 

(«)) 

(«,)) 

(«)) 
(«.)) 
(«)) 



and ffTTB into (HMlf and ffTW]) . we find 



0.132(1 
0.132(1 
0.465(1 
0.465(1 
0.374(1 
1.041(1 
1.041(1 



14.7 



J-/2 

14.7 

q-/2 

4.87 

q-/2 

4.87 
7.72 
3.42 
3.42 



J-/2 



)mH(0)2 + 
)mH(0)2 + 

)mHioy + 

,2 )"^h(0)' + 
)mHiOy - 
)mHiOY - 



gl-gl + A\K 



5 



l9l + 92 

\„2 I „2\ A 2 



A2mH(0)^ 



r2 



\g\-gl\ l^^mnW 



\9l + 



J2 



«3 I 



1^7? + gi 



59l 



q-2 ' 



(201) 



1^?? + gl 
hi + gl 



J2 

A^mj,(0)^ 
A^mH(O)^ 

r2 



Finally, using expression (11711) and taking \\u-^ = 1, these squared masses 
become 



(«)) 



(«)) 



(«)) 



(«,)) 
(«)) 



0.132m/^(0)2(l 
0.132mi^(0)2(l 
0.465mi^(0)2(l 
0.465m/^(0)2(l 
0.374mH(0)2(l 
1.041mH(0)2(l 
1.041m//(0)^(l 



14.7 


- 63.5A2 




J-2 


14.7 


-21.1A2 




q-2 


4.87 


- 3.37A2 




r2 


4.87 


- 6.69A2 




q-2 


7.72 


- 26.0A2 




q-2 


3.42 


-3.12A2 






3.42 


- 3.56A2 


r2 



0.132mjy(0)2 
0.132miy(0)2 
0.465mjy(0)2 
0.465mjy(0)2 
0.374miy(0)2 
1.04W(0)^ 
1.04W(0)2 



(202) 



Note that the squark mass squares are positive for for all ~ 



40 and 



< A < 1. The superpotential contributions to both third family up-squark 
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masses can become significant when and are simultaneously in the 
lower and upper part of their range respectively. Be that as it may, for typical 
values of these parameters such contributions are negligible. Therefore, for 
simplicity, we ignore them when giving the leading order estimate of the 
squared masses in (12021) . 

We conclude from ffTTgD. ffTT^ . ffTMj) and fl2U2]) that 

((-1,2)) = 0, ((.3)) = (1.05)- 2"^^^°) 



((e,)) = ((L,)) = z = l,2,3 (203) 



2AmH{0) 



^ yJ'l9i{tEwy + 92{tEwy 
{{H'+)) = ((^'°)) = ((^-)) = 

is a stable local minimum of the potential energy at scale tEw- 



5 The B-L/Electroweak Hierarchy 

The vacuum state (I203P spontaneously breaks both B-L and electroweak 
symmetry, and exhibits a distinct hierarchy between the two. Color and elec- 
tric charge are left unbroken. Using (I156p . we see that the S-L/electroweak 
hierarchy, expressed as the ratio of the third right-handed sneutrino and 
Higgs vacuum expectation values, is 

where the gauge parameters are computed at Iew- We have dropped the 
term proportional to 2.25(1 — A^)/tan/3^ in f ll56p since it is always much less 
than unity in our parameter regime where tanP ~ 6.32 and < < 1. 
Note that for fixed tanf] the ratio of VEVs in (12041) can be made arbitrarily 
large by fine-tuning A — > 0. Conversely, by fine-tuning A ^ 1 this ratio 
approaches 2.22 tanf3, where we used (117ip and ( 1198^ . A more natural value 
for A would lie in the middle of the range < A^ < 1. For specificity. 
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let us take A = In this case, the ratio (12041) evaluated in the region 
6.32 < tan(3 < 40 is found to be 

19.9 < < 126 . (205) 

A second measure of the i?-L/electroweak hierarchy is given by the ratio 
of the B-L vector boson mass to the mass of the Z boson. It follows from 
(IH5D,([I56D,([I80D anddHOl) that 

^ ^ (1.95)^ . (206) 

Note that if we take A ^ this mass ratio becomes arbitrarily large, whereas 

Ma 

if A — >^ 1, the upper bound in our approximation, then 'J^'^ is 1.95 tan(5. 



Mz 

- -j^ and evaluat 

in the range 6.32 < tan(3 < 40, one finds 



Again, using the more natural value A = and evaluating this mass ratio 



17.5 < < 110 . 207 

- Mz - ^ ' 

We conclude from ^^20^ and l{206\) that for typical values of A, the vac- 
exhibits a B-L/electroweak hierarchy 0/0(10) to O(IO^) over the 
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physically interesting range 6.32 < tanf3 < 40. Let us review the reasons 
for the existence and magnitude of this hierarchy. First, initial conditions 
pij) . fH2l) ■ fl63l) ■ fl6^ give emphasis to the right-handed sneutrinos by not re- 
quiring their masses be degenerate with the Lj and soft masses. This 
enables the S[ parameter fl65p not only to be non-vanishing but, in addi- 
tion, to be large enough to dominate all contributions to the RGEs with the 
exception of the gluino mass terms. This drives negative and initiates 
B-L breaking at scale m^. Second, B and M3 (hence, muitEw)) are chosen 
to satisfy constraints fll34l) . fll35l) and (11361) respectively, with < A^ < 1. 
This insures electroweak breaking for positive at a scale proportional to 
Amnio) /T. The large value assumed for T implies that the non- vanishing 
VEV is largely in the direction, allowing one to identify T, to leading 
order, with tan/3. Third, equation (11561) insures that squark/slepton squared 
masses are positive at all scales. This guarantees that the electroweak break- 
ing is substantially smaller than the B-L scale, with the 5-L/electroweak 
hierarchy proportional to tan/3/A. 
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To finish our analysis, we would like to emphasize that once the ^m~^ 
ratio is fixed, that is, once the B-L vector boson mass is measured, the 
spectrum of squarks and sleptons is completely determined. To see this, first 
recall that at Iew 

Ma^_, = V2gB^L{tEw) ((i^a)) • (208) 

It then follows from (fTSD . ffTOOj) and ffT56|) that 

mH(0) = .362 Mas_^ • (209) 

Here and henceforth we drop all terms since they are much less than 

unity in our parameter regime. It then follows from (I194p that the slepton 
masses, in the order of decreasing mass, are 

((m.,,)) = 5.26 Ma,_, , 

((m,,)) = 1.72 M^,_,, (210) 

((m^J) = ((m^J) = 1.54 M^, 

((mej) = 1.30 M^, 



Similarly, we see from f l202p that the squark masses, in descending order, are 
given by 



((m,,,,)) = ((m^dj) = .614 M^, 

{{^u,,,)) = ((mz,,,)) = .410 Ma,_, , (211) 

((m„3)) = .369 M^, 

{{t^u,)) = {{niD,)) = .219 Ma, 



^B-L ' 



Note that the sleptons masses are on the order of M4g_^ and each is heavier 
than any squark. The squark masses are lighter, being on the order of about 
20%-60% of M^^_^, with the third family left-handed up and down squarks 
being the lightest. We conclude that the radiative i?-L/electroweak hierarchy 
also leads to a computable hierarchy in the squark/slepton masses. 



Appendix A: The Numerical Solution for mH{t)'^ 

In this Appendix, we present a numerical solution of equation fllOSp . This is 
accomplished using the numerical fitting and solving packages in Wolfram's 
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Mathematica program [BB]- The details of these packages are quite complex 
and we refer the reader to the documentation provided by Mathematica for 
further discussion. 

We begin by considering equation (11081) with the property that is 
positive over the entire scaling range Iew ^ t < 0. With this additional 
condition, fllUSI) is equivalent to equation (11091) . First note, using 1^^, that 
f llOSp contains two arbitrary parameters, mj^(O)^ and |M3(0)p. The coeffi- 
cient m//(0)^ is the value of the up-Higgs soft supersymmetry breaking mass 
at t = and, hence, is a natural input parameter. However, for the reasons 
discussed in the text, it is more transparent physically to input the value of 
mjj at tEw, given in (fT^ . (fT^ and (fT^ . rather than |M3(0)p. To do this, 
one chooses T and subject to the constraints ^ 40 and < < 1. 
It then follows from fll34p and fll35p that we have completely specified the 
value of rriHitEw)^- For this choice of parameters, we solve the mjj equation 
of motion fllOSp numerically, adjusting the value of |M3(0)p until the value 
of m^fj at tEw is given by m//(0)^/T'^. The plot of this value of |M3(0)p for 
a large number of choices of T'^ is shown in Figure 1. Fitting this result with 
a smooth curve, we find that 

|M3(0)p = .0352(1 - ^)mH{Of • (212) 

This curve is also plotted in Figure 1 and closely reproduces the numerical 
data. This justifies expression (I138P used in the text. 

Having determined relation (I212p . one can now numerically solve the 
equation (llOSp for any choice of m//(0)^ and T'^. Numerical plots of 
m'jj/mH{Oy as a function of t are shown in Figure 2 for several different 
choices of T'^. For each such graph, we fit the numerical results to a smooth 
curve. We find that for any allowed value of T'^, the data is well represented 
by 

m^itr = (1 - (1 - ^)(^1^)(^))^^(0)' ' (213) 
where b is a function of T' of the form 

b(T') = 19.9(1 - ^ '^^^ ) . (214) 

Note that at t = and t = tEWi ^'h is given by m//(0)^ and mui^)'^ /T'"^ 
respectively, as it must be. This smooth curve is plotted in Figure 2 for each 
of the choices of T'^and is seen to give a close fit to the numerical data. This 
justifies equations fll39p . (ll40p used in the text. 
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Fi gUre li This plot shows a representative set of points of |A'/3(0)p/m//(0)^ (black dots) for different 
values of T'^ as well as the accuracy of the fit of equation I I212I I (blue line) to these representative points. 
For simplicity we show only show a few points but many more were used in the generation of this fit. 
The accuracy of this fit over the range of T'^ exceeds one percent at each data point (including those not 
shown). 
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(a) T - 6.32 
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(b) T' = 40 
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(c) T' = 120 

Fi gUre 2: These plots show representative points m[j{t)^ /m[j{Qi)^ (black dots) obtained by numerical 
solution of (I108I I for different values of f in our scaling range as well as the fit of equations (I213I I , (I214I I (blue 
line) to these points. We show this for three choices of T' spanning a wide range of physically interesting 
values. 
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Appendix B: The Relationship Between mi/(0)^ and mif(O)^ 

As highlighted in the discussion following equations (11551) and (11931) . squark 
and slepton masses evaluated around the VEVs {{1^3)), {{H'^)) generically de- 
pend on the two initial values, m^(0)^ and 771/^(0)^. We want to explore the 
region of parameter space leading to positive squark / slepton masses over the 
entire scaling range, thus simplifying the discussion of symmetry breaking. 
This region of parameter space can be specified in terms of a simple T'^ de- 
pendent relationship between these two initial conditions. In this Appendix, 
we present a detailed derivation of this relationship. 

Recall that the squark/slepton masses depend not only on the associated 
coefficients and their running, but on contributions from the ((^3)) and 
VEVs as well. These expectation values are only non-zero below cer- 
tain scales, thus complicating the analysis. Furthermore, both mg^ and m^^ 
contain an m'jj term and, hence, require the numerical solution described in 
Appendix A. As an explicit example of how the relationship between m^(0)^ 
and m/^(0)^ effects the sign of the scalar mass terms, let us consider ((mg^)). 
It follows from ([HHD, (^M> and (HHS]) that 

(Ka)) - 
+ 

where =F indicates that for the up field U3 and the down field of the doublet 
Q3 one uses a minus sign and plus sign respectively. Let us discuss each term 
in this equation and its value over the entire scaling range. The first term 
is m'jj derived in Appendix A and given in (12131) . (12141) . It is proportional 
to miif(O)^ and depends explicitly on the value of T'. We plotted it over 
the scaling range for three different values of T' in that Appendix. Next 
we have the integrals containing IM3P and S[. These can be evaluated for 
arbitrary t using equations (!99l) . (ll38p and (176|) . (ll00p respectively. Note that 
the first is T'^ and m//(0)^ dependent, while the second is m;^(0)^ dependent. 
Figure [3] plots their values over the scaling range, where in graph (a) we have 
taken T' = 40 for specificity. Next we have the sneutrino VEV term. It is 
given for any scale t by minimizing the purely z/3 part of potential (1771) using 
equations f|T9l) . fl57j) and fl76|) respectively and is plotted in Figure H] over the 
scaling range. Note that this term depend on m^(0)^. Lastly, we have the 
Higgs VEV term. This turns out to give a very small contribution, so we 
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Figure 4: in this graph, ((i/3))2/m,.(0)2 is plotted over the entire scaling range. 

omit it henceforth. 

Note that each of the above terms is proportional to either m^lOY or 
'^//(O)^. Furthermore, three of the terms are everywhere positive whereas 
the two remaining terms, which are proportional to m,y(0)^, are everywhere 
negative. It is easy to demonstrate that for a sufficiently large ratio of mj,(0)^ 
to m/f(0)^, {{fng,^)) could be negative somewhere in the scaling range. To 
investigate this, we first define 

m^(0)2 = D^mniOY ■ (216) 

This allows us to write each term as a function of one of these initial masses, 
which we choose here to be m^(0)^, and their ratio D^. Inserting the renor- 
malization group expression for the scaling of each term into fl215p . we find 
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that 



(Ka)) ^ o (1 - (1 - :^)(^-— )(r3l))- 



1)0.0352 



11.5) m,(0) 

T-/2 ' n2 



^ - 1)149 m^(0)2 ^ ' ^ 



m.(0)^ - ^(1 - (1 - ^^^)-W)149m.(0)^ (217) 

m ( -2 



4(l-^)H^7(0)(64-i) 
(7(0)%t\5(i-2li) 



87r2 



149m^(0)^ 



We have preserved the original ordering of terms for ease of reference. Recall 
that the coefficient h is defined in f l214p . Note that we can factor out m,^(0)^ 
and, hence, the properties of this equation are controlled by the value of 
and T'^. Choose a fixed value for T'^. If we initially assume that <C 1, 
that is, m,^(0)^ ^ mj|/(0)^, then {(juq^) is everywhere positive. However, 
as we increase the value of ((mg^)) becomes smaller, eventually touching 
zero at some point t in the scaling regime. For still larger D"^ the value of 
((mg^)) at this point, and in a range around it, is negative, see Figure 
Record the value of D"^ at which {{rriQ^)) just touches zero and repeat this 
for a set of choices of T'^. This is plotted in Figured We find that a smooth 
curve fitting this data is given by 

2 25 

= 0.864(1 - ^) (218) 
or, equivalently, using f l216p that 



2 25 

m,(0)2 = 0.864(1 - -^)mH{OY . (219) 
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<< m^Q3 (t) >> 
IIlH (0)" 




t 



Fi gUre 5; in this graph, wc demonstrate the dependence of the running of ((mQ^)) on D^. We plot 
({rrtQ^)) /'iTiH {0)'^ for D=0.75 (green line), D=0.85 (blue line), and D=0.95 (red line). It is apparent from 
this plot that as D increases, {{rriQ )) goes negative. For this plot, we took T'^ = 40. 
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Fig U.r6 6i This plot shows representative points (black dots) for different values of T"'"^ and the 

fit of equation I I218I I (blue line) to these points. The accuracy of this fit over the range of T'^ exceeds one 
percent at each data point. 
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This equation specifies the largest value of m^^O)^ relative to m//(0)^ 
for which ((^tIq^)) is everywhere non-negative. An investigation of all other 
other squark and slepton mass squares evaluated around the B-L and Higgs 
VEVs shows that they are everywhere positive whenever {{mQ,J) is every- 
where non-negative. That is, for m^(0)^ and mji^(O)^ satisfying relation fl219p . 
all squark/slepton masses are everywhere positive, as desired. Furthermore, 
this is the largest value of m,y(0)^ relative to 771^(0)^ for which this is the 
case. This is chosen so that the B-L /EW hierarchy, which is proportional to 
my(0)/miy(0), is as large as possible within this context. This justifies our 
use of equation (11561) in the text. 



Appendix C: Comparison to the Standard Formalism 
and the Higgsino Mass 

In this Appendix, we compare our analysis of the Higgs section with the 
standard MSSM results given, for example, in "A Supersymmetric Primer" 
by Stephen Martin [6lj. Restoring the fi paramter, that is, not assuming it 
is necessarily sub-leading, the complete potential for the neutral Higgs fields 
becomes 

V = (|/i|2 + m^)|/JY + (l/^l' + "^|)I^T-5(^°^° + M 
8 5 

Note that this potential is written in terms of the original Higgs fields, not the 
mass eigenstates at the origin we used in our analysis. The two constraints 
that must be satisfied to have a stable, non-vanishing vacuum of (12201) are 
presented as equations (7.3) and (7.4) in the Supersymmetric Primer. The 
two minimization equations are given in (7.10) and (7.11) of that review. 
Note that each of these equations should be evaluated at the electroweak 
scale tEw- We now show, to the order in T^^ we are working to, that our 
vacuum satisfies all four of these equations. 

We begin by examining the constraint (7.3) in [Hlj. This equation, which 
ensures that the potential is bounded from below, and is given by 

2fi < 2|/i|2 + + m| . (221) 

Using (11361) . this can be rewritten in the form 

:^ = l/^r + ^^ + ^. (222) 
T — mjj 2 
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It then follows from fll33p . fll34p . 01351) and 01361) that this inequality is satis- 
fied by our vacuum over the entire range T ~ 6.32 and < < 1 for any 
value of the parameter /i. Next, consider the constraint equation (7.4) in [61] 
given by 

B'<i\fi\' + mj,){\fi\' + ml). (223) 

This ensures that the origin will not be a stable minimum of the potential. 
Using ffT33D.ffT3D. fn:35|) and IHM) this becomes 

1 > (1 + ^)(1 + ^77^)7-^zz^^ • (224) 

For the time being, take fi to be sub-dominant to the Higgs mass parameters 
and ignore it, as we did in the text. Later in this Appendix we will keep 
it in our analysis and observe the consequences. We are left with a simple 
inequality 

(1 - ^-j^r > 1 - (225) 

We want to examine if this holds over the entire range < A^ < 1 and 
T ~ 6.32 of our vacuum solution. Note that 0225P is trivially satisfied as 
A^ — i> 1. However, in the A^ — > limit this inequality clearly does not hold. 
To find what bound is set by the above equation, solve 0225p for A^ keeping 
the relevant root. We find that the inequality will be satisfied for 

A^ > + ^ + 0(ij) (226) 

for any value of T. To the order we are working in this paper, is in fact 
the lower bound and inequality 0225P is satisfied over the whole range of 
A^. However, to next order in T~^, a constraint appears. This is due to 
the fact that we only took the first order approximation in equation 013ip . 
Had we kept higher orders, then relation 01371) would be more complicated 
and, hence, the lower bound on A^ different. It is straightforward to verify 
that, to higher order in T~^, the lower bound would be identical to (I226p . 
We conclude that, to the order we are working, constraint (7.4) in [6l] is 
satisfied over the entire range T ~ 6.32 and < A^ < 1 of our vacuum. 

We now show that identities (7.10) and (7.11) in |Mj, which are derived 
from the minimization conditions of potential O220p . are also satisfied by our 
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vacuum. First consider (7.10) given by 



IB 

sin2(3 = 2 • (227) 





+ m\ 


+ 2|/iP 


as 








2 




C(l- 


-(1- 






B 





i^^P = 4TTm ' (228) 

where 

= [ml + ml + 2\fi\^ ) ■ ^^^^^ 
Expressing C in terms of T using f ll36p . we find that f l228p becomes 

tanP = r(l + ^ + . (230) 

Hence, to leading order, T = tan(3. This is identical to expression ( 1183^ in 
our vacuum solution. Note that this is true for any value of the parameter 
fi. It also follows from f l230p that to higher order in T~^, the relationship of 
T to tanf3 becomes non-linear. Finally, consider the minimization equation 
(7.11) in [6l]. For our purposes, this is most conveniently re-expressed for 
large tanf3 in (7.12) as 

,4 = -2i^h + \,?) + ^(,4 - <) + O(j^) . (231) 

It is clear from fll33p . fll34p . (11801) ■ fll83p that as long as fi is chosen to be sub- 
dominant to niH, our vacuum explicity solves fl23ip . We conclude that, to 
the order we are working, minimization equations (7.10) and (7.11) in [61] 
are satisfied over the entire range r ~ 6.32 and < < 1 of our vacuum. 

An assumption in this paper was that /i, while non-vanishing, was sub- 
leading to the Higgs mass parameters at tEw and, hence, ignorable in the 
vacuum solution to leading order. Can one quantify how much this assump- 
tion restricts the size of /i and, hence, physical quantities such as the Higgsino 
mass? To do this, note from this assumption and fll33p . fll34p and fll35p that 

1 - 

< ~ __m^(0)2 < m| ~ mniO) (232) 
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in our vacuum. Therefore, a reasonable parameterization of /i which is sub- 
leading to the Higgs masses is 

1^1^ "(^-^^'. (233) 



(234) 



ml T2 

where a is a constant of (9(1). It then follows from (I134p that 

l/ip _ a(l- A^)^ 

Putting this into inequality (12241) . one finds the condition 

A' > (235) 

to the first non-vanishing order in T^^. To see how this limits the size of n, 
calculate the ratio of to Mz using (11801) and (12331) . We find 

"I - (236) 



Mz V 2 A T 
Applying constraint (I235p . it follows that 



Mz V2(" + 2)' T2 



which, for the range T ~ 6.32 and taking, for example, a = 1, gives 

l^il < -^Mz = 37.5GeV . (238) 
v6 

We conclude that, even subject to our assumption that it be sub-leading 
to the Higgs mass parameters at tEw, the /i parameter and, hence, such 
quantities as the Higgsino mass, can be relatively large, with an upper bound 
approaching 40% of the Z-mass. 
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